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WHAT SHAPES THE STRUCTURE OF MOLECULAR CLOUDS: TURBULENCE OR GRAVITY? 
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ABSTRACT 

We revisit the origin of Larson's scaling relations, which describe the structure and kinematics of molecular 
clouds, based on recent observations and large-scale simulations of supersonic turbulence. Using dimensional 
analysis, we first show that both linewidth-size and mass-size correlations observed on scales 0.1-50 pc can 
be explained by a simple conceptual theory of compressible turbulence without resorting to the often assumed 
virial equilibrium or detailed energy balance condition. The scaling laws can be consistently interpreted as a 
signature of supersonic turbulence with no need to invoke gravity. We then show how self-similarity of structure 
established by the turbulence breaks in star-forming clouds through development of gravitational instabilities 
in the vicinity of the sonic scale, 4 ~ 0.1 pc, leading to the formation of prestellar cores. 
Subject headings: stars: formation — ISM: structure — turbulence — methods: numerical 



1. INTRODUCTION 

iLarsonl (Qj981) established that for many molecular clouds 
(MCs) their internal velocity dispersion, <j„, is well corre- 
lated with the cloud size, L, and mass, m. Since the power- 
law form of the correlation, a u o c L 038 , and the power in dex, 
0.38 ~ 1 /3, were similar to the Kolmogoro^ (1 1 94 1 albT) law 
of incompressible turbulence (K41), he suggested that the ob- 
served nonthermal linewidths may originate from a "common 
hierarchy of interstellar turbulent motions." Larson also no- 
ticed that the clouds appear mostly gravitationally bound and 
in approximate virial equilibrium, as there was a close posi- 
tive correlation between the velocity dispersion and mass of 
the clouds, a u oc m 0,20 , but suggested that these structures 
"cannot have formed by simple gravitational collapse" and 
should be at least partly created by supersonic turbulence. 
This seminal paper preconceived many important ideas in the 
field and strongly influenced its development for the past 30 
y ears. 

ISolomon et alj d 1987ft confirmed Larson's study using ob- 
servations of I2 CO emission with improved sensitivity for a 
large uniform sample of 273 nearby clouds. Their linewidth- 
size relation, tj u = 1 .0 ± 0. 15 ,0 5±0 05 km s" 1 , however, had a 
substantially steeper slope than Larson's, reminiscent of that 
for clouds in virial equilibrium^ 

a u = (nGZ) l/2 R 1 / 2 , (1) 

since the ISolomon et al.l d 19871) clouds had approximately 
constant molecular gas surface density, E, independent of 
their radius. The surface density-size relation, also known as 
the third Larson's law, can be derived eliminating er„ from his 
first two relations: E oc mL~ 2 oc ^°- 38 / - 2 °- 2 (x ZT° 1 Assuming 
that E = const for all clouds, Solomon et al. ( 1987) evaluated 

1 We deliberately keep the original notation used by different authors for 
the cloud size (e.g., the size parameter in parsecs, S = Dtan(^/<r/(76); the max- 
imum projected linear extent, L; the radius, R = ^A/ir, defined for a circle 
with area, A, equivalent to that of cloud) to emphasize ambiguity and large 
systematic errors in the cloud size and mass estimates due to possible lin e-of- 
sight confusion, ad hoc cloud boundary definitions i Heyer et al. 2009), and 
various X-factors involved in conversion of a tracer surface brightness into 
the H2 column density. 



the "X-factor" to convert the luminosity in 12 CO (1-0) line to 
the MC mass. The new power index value ^0.5 ruled out Lar- 
son's hypothesis that the correlation reflects the Kolmogorov 
law a u oc L 1 / 3 . In the absence of robust pr edictions for the ve - 
locity scaling in supersonic turbulence (cf. lPassot et aljl9"88l) . 
simple virial equilibrium-based interpretation of linewidth- 
size relation appealed to many in the 1980s. 

Since then views on this subject remain polarized, while 
it is still uncertain whether gravity or turbulence (or mag- 
netic fields) govern the dynamic s of MCs. For instance, 
iBalles teros-Paredes et al l d201 latfbh argue that MCs are in 
a stat e of "hierarchical a nd chaotic gravitational collapse," 
while iDobbs et al.l d201 ll) believe that GMCs are "predomi- 
n antly gravitation a lly un bound objects." 

Heyer& Brunt (2004) found that the scaling of velocity 
structure functions (SFs) of 27 GMCs is near invariant, 



S l (u,£)={\Su e \) = u £ > 



0.56±0.02 



(2) 



for cloudy structures of size £ 6 [0.03,50] pcQ In the mean 
time, numerical simulations of supersonic isothermal turbu- 
lence returned very similar inertial range scaling exponents 
for the first-order velocity SFs (Si(u, £) oc £^ with Ci = 0.53 ± 
0.02 and 0.55 ± 0.04 for longitudinal and transverse SFs, re- 
spectively, see Fig. Q] and Kritsuk et al. 2007a). The rms 
sonic Mach number, M s = 6, used in these numerical experi- 
ments is characteristic of MC size scale £ w 2 pc, i.e. right 
in the middle of the observed scaling range. The simula- 
tions indicated that the velocity scaling in supersonic regimes 
deviates strongly from Kolmogorov's predi ctions for fluid 
turbulence. This result removes one of the ISolomon et alJ 
(119871) arguments against Larson's hypothesis of turbulent 
origin of the linewidth-size relatio n. Indeed, the scaling ex- 
ponent of 0.50 ± 0.05 measured by ISolomon et all (119871) for 
whole clouds and a more recent a nd precise measurement 
0.56 ± 0.02 bv lHever & BruntJ (120041) that includes cloud sub- 
structure both fall right within the range of expected values for 

2 The lengths entering this relation are the characteristic scales of the PCA 
eigenm odes, therefore they ma y differ from the cloud sizes defined in other 
ways iMcKee & Ostriker 2007). 
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Figure 1. Scaling of the first-order transverse (green) and longitudinal (red) 
velocity SFs in a sim ulation of isothermal supersonic turbulence with M s = 6 
I Kritsuk et al. 2007a). A is the grid spacing. 

supersonic isothe rmal turbulence at r elevant Mach numbers 
More recently, iHever et all (|2009) used observations of a 
lower opacity tracer, 13 CO, in a sample of 162 MCs with 
improved angular and spectral resolution to reveal system- 
atic variations of the scaling coefficient, Uq, in Equation (O 
with I and S. Motivated by the concept of clouds in self- 
gravitating equilibrium, they introduced a new scaling coeffi- 
cient u' = (\Su£\)i~^ 2 oc S°' 5 . This correlation would indicate 
a departure from "universality" for the velocity SF scaling © 
and compliance with the virial equilibrium condition ((T). 

An alternative formulation of the original Larson's third 
law, mrxL 19 , implied a hierarchical density structure in MCs. 
Such concept was proposed 60 years ago by von Hoerner 
to describe a complicated statistical mixture of shock waves 
in highly compressible inter stellar turbulence (Ivon Hoerneri 
1195 lb Ivon W eizsackerl 119511) . Von Hoerner pictured den- 
sity fluctuations as a hierarchy of interstellar clouds, anal- 
ogous to eddies in incompressible turbulence. MC obser- 
vations indeed reveal a pervasive fractal structure in the in- 
terstellar gas that is interpre t ed as a signature of turbulence 
(iFalgarone & Phillips! 11991b lElmegreen & Falgarond 11996b 
Roman-Duval et al. 2010). The most r ecent result for a sam- 
ple of 580 MCs, which includes the ISolomon etail (119871) 
clouds, shows a very tight correlation between cloud radii and 
masses, 

m(R) = (228 ± 1 8M )fl 2 ' 36±004 , (3) 

for R e [0.2,50] pc (iRoman-Duval et al.ll2Q10h . The power- 
law exponent in this relation is simply the mass dimension 
of the clouds, d m w 2.36, which corresponds to a "spongy" 
medium organized by turbulenc e into a multiscale pattern 
of clustered corrugated shocks (iKritsuk et al.l [2006). Di- 
rect measurements of d m in high-resolution three-dimensional 
simulations of supersonic turbulence give the inertial sub- 
range values in excellent agreement with observations: d m = 
2.39 ± 0.01 (1024 3 grid cells. lKritsuketal.l2007al) and 2.28 ± 
0.01 (2048 3 cells, F ig. EJ, depending on details of forcing 
(IKritsuk et al.l |2010). Note, that d m pa 2.36 implies £ oc 
mL~ 2 oc L 36 , thus the observed mass-size correlation does 
not support the idea of a universal mass surface density of 
MCs. Meanwhile, positive correlation of £ with £ removes 

3 Note that the distinction between scaling "inside clouds" and "between 
clouds" often used in observational literature is specious because clouds are 
not isolated entities on scales of interest (e.g., Henriksen & Turner 1984). 



Figure 2. Compensated scaling of the mass , me, with size I in a 2048 3 simu- 
lation of isothermal turbulence with M s = 6 ( Kritsuk et al. 2009). Dissipation- 
scale stru ctures are shocks with d m ft: 2. Inertial-range structures have 
d m w 2.3 iKritsuk et alJ2007al 1200911 . 

theoretical objections against the third Larson's law outlined 
above. 

With d m pa 2.36, the power-law index in the linewidth-size 
relation compatible with the virial equilibrium condition ([T), 
Ci vh- = (d m - 1)/2 ps 0.68, is still reasonably close to the scal- 
ing exponent (i ~ 0.56 in Eq. (f2|, even if one assumes uq = 
const (see §3 below). Thus we cannot immediately exclude 
the possibility of virial equilibrium (or kinetic/g ravitational 
energy equipartition, see Balle steros-Paredesl2006l) across the 
scale range of up to three decades based on the available ob- 
servations alone. At the same time, we have seen that the 
two classes of observed correlations (e.g., linewidth-size and 
mass-si ze) are readily reprod uced in simulations without self- 
gravity (Kritsuk et al. 2011b). In this sense, non-gravitating 
turbulence is self-sufficient at explaining the observations. 
Hence, following Occam's razor, it is not necessary to in- 
voke gravity. Meanwhile, in turbulence simulations with self- 
gravity, the velocity power spectra do not show any signa- 
ture of ongoing core formation, while the density and column 
densit y statistics bear a s trong gravitational signature on all 
scales (Co llins et alj20 12). What is the nature of this apparent 
"conspiracy" between turbulence and gravity in MCs? Why 
do structures in MCs appear gravitationally bound when they 
might not be? 

In this Letter we use a simple conceptual theory of super- 
sonic isothermal turbulence to show that scaling exponents in 
the linewidth-size and mass-size relations are connected, i.e. 
one can be derived from the other. In §2 we briefly introduce 
the concept of compressible cascade and discuss potentially 
universal relations. In §3 we derive the surface density-size 
and m,')-S relations in several different ways and demonstrate 
consistency with observations and numerical models. §4 deals 
with the effects of self-gravity on small scales in star-forming 
clouds, and discusses the origin of the observed mass-size re- 
lation for prestellar cores. Finally, in §5 we formulate our con- 
clusions and emphasize the statistical nature of the observed 
scaling relations. 

2. WHAT'S UNIVERSAL AND WHAT'S NOT 

In turbulence research, universality is usually defined as in- 
dependence on the partic ular mechan ism by which the turbu- 
lence is generated (e.g., Frisch 1995). Following this conven- 
tion, the non-universal nature of scaling relation (O can be 
readily understood. Indeed, any scaling law for compressible 
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Figure 3. Absolute scaling exponents for transverse SFs of velocity 
(S p (u,£) <x ISp, circles) and density-weighted velocity (S p (v,t) oc 

t^p , crosses) from Kritsuketal. 1 2007a). Solid lines show Q, predicted by 
th e K41 theory (red), Bu rgers' mode l (blue), and inte rmittency models due 
to She & Leveque 11994, green) and Boldyrev 12002, magenta). Note that 
Burgers' model predicts Ct = 1 (cf.[McKee & Ostriker 2 0(3% . 

turbulence formulated in terms of the velocity alone would 
depend on the Mach number. At M s < 1, density fluctua- 
tions are relatively small and turbulence re mains very similar 
to the inco mpressible case with Ci ~ 1/3 dPorter et al.ll200l 
iBenzi et aLll2008l) . The scaling steepe ns at higher Mach num - 
bers reaching Ci ~ 0.54 at M s sa 6 (iKritsuk et al.1 l2007allbt 
iPan & Scannapiecoll201 lb . This transition is also accompa- 
nied by a change in the dimensionality of the most singular 
dissipative structures from 1 (vorte x filaments) to 2 (shocks, 
iBoldvrevllIOOllPadoan et aHl2004l) . 

In order to get a universal scaling relation, one would usu- 
ally have to consider ideal quadratic invariants or, more pre- 
cisely, rugged invariants. In compressible isothermal flows, 
however, the situation in more complex. There are two invari- 
ants: the mean total energy density E = (pu 2 /2+c 2 pln(p/ po)) 
(po is the mean density) and the mean helicity. Since pu 2 /2 is 
positive semidefinite, its contribution to E dominates at high 
Mach numbers, while the second term is always subdominant. 
Hence, in the case of MC turbulence, the kinetic energy den- 
sity is probably the best possible choice (if any), as far as 
universality is concerned. In MCs, on scales above the sonic 
scale, I ^ £ s ~ 0.1 pcQ turbulence is highly compressible and 
any quantity with universal scaling must depend on both mass 
density and velocity. 

An illustration is given in Figure [3] where filled circles 
show the absolute scaling exponents, ( p , of the velocity SFs 
of order p £ [0.5,5.5] at M s = 6. These exponents would be 
close to the K41 prediction, ( p = p/3, at M s < 1, if we as- 
sumed isotropy, homogeneity, and ignored intermittency cor- 
rections. At M s J> 3, C3 shows an excess over unity, which 
systematically increases with the Mach number, indicating a 
non-universal behavior. A better candidate for universal scal- 
ing is the density-we ighted velocity, v = p x J^u since £3 « 1 
at all Mach numbers (Rri tsuk et d]l2007allbh . We will use the 
third-order moment of Sv as its linear scaling implies approx- 
imately constant kinetic energy flux within the inertial range, 

S 3 (v,£)=(\Sv ( \ 3 ) = (e)£, (4) 
with some minor intermittency correction (Kritsuketal] 

4 The sonic scale £ s is defined by the condition = c s . 
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Figure 4. Variation of the scalin g coefficient u' n = a v R' l l 2 with mass surface 
density E based on data from Heyer et al. 1 2009). Solid lines with slopes 
0.34 ±0.04 and 0. 32 ± 0.03 show least square fits to Al and A2 subsets from 
Heyer et al. ( 2009), respectively. Dashed line shows the correlation expected 
for clouds in virial equilibrium. 

I2007at lAlui3l20Tlt iGaltier & Banerieell20Tlb . Here, (e) is 
the density-weighted energy transfer rate independent of I. 
Lower and higher than p = 3 moments of 5v would most prob- 
ably require sub stantial intermittenc y corrections that are hard 
to predict dLandau & Lifshitzll 19441) . 

3. EXPOSING MOLECULAR CLOUD CONSPIRACY 

In the following, we exploit the "universal" third-order rela- 
tion to derive several secondary scaling laws involving S^, 
assuming that the effective kinetic energy flux across the hi- 
erarchy of scales is approximately constant within the inertial 
sub-range (and neglecting gravity). 

Dimensionally, the constant spectral energy flux condition, 

p e (Sue) 3 t x oc E £ (Su e f t 2 oc Hd K '~ 2 ~ const, (5) 
implies 

S £ oc £ 2_3Cl . (6) 

Substituting Ci = 0.56 ± 0.02, as measured bv lHever & Brunj 
(2004), we get S £ oc £ 32 ±o.06 We can also rely on me f racta i 
properties of the density distribution to evaluate the scaling of 
T,£ with t. oc pil (xmi/l 2 oc £ dm ~ 2 , which in turn implies 
(xj 0-36±0-04 f or d m = 2.36 ± 0.04 from lRoman-Duval et alJ 
(2010). Note that both independent estimates for the scaling 
of with £ agree with each other within one sigma. New 
observational data, thus, indicate that mass surface density of 
MCs indeed positively correlates with their size with a scal- 
ing exponent ^1/3. This result is consistent with both ob- 
served velocity scaling and observed self-similar structure of 
the mass distribution in MCs. 

Let us now examine data sets Al and A2 presented in 
iHever etaT] d2009t) for a possible correlation of u' = a u R~ x l 2 
with E = m/nR 2 . Figure [4] shows formal least-square fits for 
the two data sets with slopes 0.34 ± 0.04 and 0.32 ± 0.03, re- 
spectively. Note that both correlations are not as steep as the 
virial equilibrium condition (Q~|) would imply. When the two 
data sets are plotted together, however, the apparent shift be- 
tween Al and A2 points (caused by different cloud bound- 
ary definitions in Al and A2) creates an impression of virial 
equilibrium condition (dash ed line in Fi g. [4]i being satisfied, 
although with an offset that lHever et alJ d2009t) interpret as a 
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consequence of LTE-based cloud mass underestimating real 
masses of the sampled clouds. Each of the two data sets, 
however, suggest scaling with a slope around 1 /3 with larger 
clouds of higher mass surface density being closer to virial 
equilibrium than smaller struc t ures. The same tendency can 
be traced i n the iBolatto et ap (120081) sample of extragalac- 
tic GMCs dHever et al.ll2009l Fig . 8). A similar trend is re- 
covered by Goodma n et al.l (120091) in the LI 448 cloud, where 
a fraction of self-gravitating material obtained from dendro- 
gram analysis shows a clear dependence on scale. While most 
of the emission from the L1448 region is contained in large- 
scale bound structures, only a low fraction of small-scale ob- 
jects appear self-gravitating. 

Let us check a different hypothesis, namely whether the 
observed scaling cr u R~ ~ l l 2 oc S 1 / 3 is compatible with the tur- 
bulent cascade phenomenology and with the observed fractal 
structure of MCs. The constant spectral energy flux condi- 
tion, pi (Sue) 3 £~ l rs const, can be recast in terms of T,£ oc pgl 
assuming Suet 1 ! 2 oc with a w 1/3 and pe oc £ dm ~ 3 , 

Pi (5u e ) 3 r 1 oc piEi&Pt 1 oc e 2(dm - 3)+3/2 » const. (7) 

This condition then simply reads as 2(d m - 3) + 3/2 w or 
d m ~ 2.25, which is close to the observed fractal dimension 
of MCs. 

As we have shown above, the observed correlation of the 
scaling coefficient u' with the coarse-grained mass surface 
density of MCs is consistent with a purely turbulent nature 
of their hierarchical structure and does not require any addi- 
tional assumptions concerning virial equilibrium. The origin 
of this correlation is rooted in highly compressible nature of 
the turbulence that implies density dependence of the lhs of 
equation Let us rewrite dU for the first-order SF of the 

density-weighted velocity: (|cH^|) ~ {e 1 / 3 )? 1 / 3 . Due to inter- 
mittency, the mean cubic root of the dissipation rate is weakly 

scale-dependent, (e]/ 3 ) oc £ Tl / 3 , and thus (|<5t>g|) oc ft 1 , where 
£i = I/3 + T1/3 and n/3 is the intermittency correction for the 
dissipation rate. Using dimensio nal arguments, on e can ex- 
press the scaling coefficient in the lHever et al. (2009) relation, 

Su e r l/2 oc pf' 3 t l l^ oc s; I/3 ^/ 6+r '/3. (8) 

Since, as we have shown above, £^ oc I 1 ' 3 , one gets 



(9) 



Numerical experiments gi ve Tin ~ 0.055 for the density- 
weighted dissipation rate dPan et al.l 120091) . This value im- 
plies s caling, Sue£~ 1 ^ 2 oc E 33 , consistent with the lHever et al.l 
(120091) data. 

4. A PLACE FOR GRAVITY 

So far, we limited the discussion of Larson's linewidth- 
size and mass-size relations to scales above the sonic scale 
£ s ~ 0.1 pc. Theoretically, the linewidth-size scaling index is 
expected to approach (1 « 1/3 at £ < 4 in MC sub-structure s 
not affected by s elf-g ravity (see §2 and Kritsuk et al. 2007a). 
Falgar one et al.l (120091) explored the linewidth-size relation 
using a large sample of 12 CO structures with £ £ [10~ 3 , 10 2 ]pc. 
These data approximately follow a power law Sue oc £ 1 ' 2 for 
£ > lpc. Although the scatter substantially increases be- 
low lpc, a slope of 1/3 "is not inconsistent with the data." 
12 CO and I3 CO observations of translucent clouds indicate 
that small-scale structures down to a few hundred AU are 



possibly intrinsically linked to the forma tion process of MCs 
(Fal garone et al.lll998l:lHeithausenll2004l) . 

The observed mass-size scaling index, d m « 2.36, is ex- 
pected to remain constant for non-self-gravitating structures 
down to £ v ~ 307/, which is about a few hundred AU, assum- 
ing the Kolmogorov scale 77 ~ 10 cm (lKritsuketal.ll2011d) . 
This trend is traced down to ~ O.Olpc with recent Herschel 
detecti on of ~ 300 unbou nd starless cores in the Polaris Flare 
region ( lAndre et al.l l2.Q10l). For scales below £ v ~ 200 AU, in 
the turbulence dissipation range, numerical experiments pre- 
dict convergence to d m ~ 2 due to shocks, see Fig. [2] 

In star-forming clouds, the presence of strongly self- 
gravitating clumps of high mass surface density breaks self- 
similarity imposed by turbulence. One observational signa- 
ture of gravity is the build-up of a high-end power-law tail 
in the column density PDF associated with fi lamentary struc- 
tures harboring prestellar cores and YSOs (IKainulainen et al.l 
12009b lAndre et all 1201 11) . The power index of the tail, p = 
— 2/(ti— 1), is determined by the density profile, p oc r~", of 
a stable attractive self-similar collapse solution ap propriate to 
the sp ecific conditions in the turbulent cloud (Krits uket al.l 
1201 lal) . In numerical simulations with self-gravity, we inde- 
pendently measured p » 2.5 and n w 1.8 in agreement with 
the theoretical prediction. This implies d m = 3 - n w 1 .2 for 
the mass-size relation on scales below ^0.1 pc. Mapping 
of the active star-forming Aquila field with Herschel gave 
p = 2.7±0.1 and d m = 1.13 ±0.07 for a sample of 541 star- 
less cores with deconvolved FWHM siz e £ £ [0.01,0.1] pc 
(iKonyves et all 120101: lAndre et all 120111) . Using the above 
formalism, we get d m = 3-n = 2 + 2/p rs 1.26 in reasonable 
agreement with the direct measurement. Overall, the expected 
mass dimension at scales where self-gravity becomes dom- 
inant should fall between d m = 1 (Larson-Penston solution, 
n = 2) and d m = 9/7 ~ 1 .29 (pressure-free collapse solution), 
see lKritsuk et all (1201 lab . 

The characteristic scale where gravity takes control over 
from turbulence can be predicted using the linewidth-size and 
mass-size relations discussed in previous §§. Indeed, in a 
turbulent isothermal gas, the coarse-grained Jeans mass is a 
function of scale £ 



J 3 -1/2 

m\ oc aePe °c 



p~ l/2 oc P~ d ^ 2 oc £ 032 if ^4 

5v 3 pf 2 OC £l+3(3-^)/2 K £1.96 if£ > 



where a 2 = Su 2 + c 2 dChandra sekhar 1951) and we assumed 
that erf w c 2 at £ < £ s , while <rf « Suj at £ ^ £ s . The di- 
mensionless stability parameter, pe = nie/m J e , shows a strong 
break in slope at the sonic scale: 

f t^-^oct 2 - 04 if£<£, 

and a rather mild growth above £ s . Since both pe and the free- 
fall time, 



tf =[3ir/(32Gpe)] 1/2 (x £ (3 ~ d ' 



V 2 



M.32 



(10) 



correlate positively with £, a bottom-up nonlinear develop- 
ment of Jeans instability is most likely at I J> £j, where pe s = ^- 
Note that both pe and tf grow approximately linearly (i.e. rel- 
atively weakly) with at £ > £ s , while below the sonic scale 
peotY*^. This means that the instability, if at all present, shuts 
off rather quickly below £ s , i.e. £j ~ £ s . The formation of 
prestellar cores would be possible only in sufficiently over- 
dense regions on scales around £ s ^0.1 pc. The sonic scale, 
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thus, sets the characteristic mass of the core mass function, 
m£ s , an d the threshold mass surface density for star for mation, 
£4 fcf. lKrumholz & McKedl2005t lAndre etai]|20Toh . 

5. CONCLUSIONS AND FINAL REMARKS 

We have shown that, with current observational data for 
large samples of Galactic MCs, a modern version of Larson's 
relations on scales 0.1-50 pc can be interpreted as an empir- 
ical signature of supersonic turbulence fed by the large-scale 
kinetic energy injection. Our interpretation is based on the 
phenomenology of highly compressible turbulence and sup- 
ported by high-resolution numerical simulations. 

Gravity can nevertheless help accumulate the largest molec- 
ular structures ( > 50 pc) that appear gravitationally bound. 
High-resolution simulations of cloud formation in the general 
ionized/atomic/molecular turbulent ISM context are needed 
to demonstrate that molecular structures identified as bound 
in position-position-velocity space are indeed genuine three- 
dimensional objects. On small scales, in low-density translu- 
cent clouds, self-similarity of turbulence can be preserved 
down to ~ 10~ 3 pc, where dissipation starts to become im- 
portant. In contrast, in overdense regions, the formation of 
prestellar cores breaks the turbulence-induced scaling and 
self-gravity assumes control over the slope of the mass-size 
relation. We show that the transition from turbulence- to 
gravity-dominated regime in this case occurs close the sonic 
scale 4 ~ 0.1 pc, where structures turn gravitationally unsta- 
ble first, leading to the formation of prestellar cores. 

Our approach is essentially based on dimensional analysis 
and the results are valid in a statistical sense. This means 
that the scaling relations we discuss hold for averages taken 
over a large number of statistically independent realizations of 
a turbulent flow. Relations obtained for individual MCs and 
their internal substructure can show substantial statistical vari- 
ations around the mean. The scaling exponents we discuss or 
derive are usually accurate within rj (5— 10)%, while scaling 
coefficients bear substantial systematic errors. Homogeneous 
multiscale sampling of a large number of MCs and their inter- 
nal structure (including both kinematics and column density 
mapping) with CCAT, SOFIA and ALMA will help to detail 
the emerging picture discussed above. 

We thank Philippe Andre for providing the scaling expo- 
nents based on Herschel mapping of the Aquila field. This 
research was supported in part by NSF grants AST-0808184, 
AST-0908740, AST-1 109570, and by TeraGrid allocation 
MCA07S014. We utilized computing resources provided by 
SDSCandNICS. 
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